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ABSTRACT 

The thermal and suprathermal processes involving solar wind electrons 
are discussed from a theoretical point of view. A model for the electron 
distribution function, f e , based on the solutions of the Boltzmann equation 
in Krook's approximation is outlined: the angular and energy dependences of 

f e for various distances from the sun between the coronal base and the earth 
are presented with the express purpose of obtaining the radial profile of 
the heat flux generated by the suprathermal s. The calculations take into 
account the effects of the morphology of the interplanetary magnetic field, 
the electrostatic polarization potential and coulomb scattering, but assume 
that, as far as the suprathermal s are concerned, the effects of wave- 
particle interactions are negligible. This surmised knowledge of the heat 
flux allows us to solve the basic fluid- dynamical equations for the solar 
wind plasma as a whole along any given magnetic field line. The results, 
although containing a variety of idealizing assumptions, are in surprisingly 
good agreement with observations. In particular, the predicted density and 
temperature profiles for positive ions exhibit the characteristics of recent 
measurements, both in the corona (above the coronal holes) and between 0.3 
and 1 A.U. (the Helios region). It is concluded that the suprathermal 
electrons play an important role in the formation and the dynamics of the 
solar wind. 


1. Introduction 


The discovery of Coronal Holes confirmed a belief of long standing that 
the stream structure of the solar wind in the steady state is closely 
related to the morphology of the magnetic field lines in the corona. The 
open field lines in the regions above the Coronal Holes serve as a nozzle 
for the coronal plasma to flow into the far regions of the Heliosphere. The 
high-velocity streams observed by various spacecraft can be thought of as 
originating in the central portion of a Coronal Hole above the coronal base 
(say, between 1.03 r 0 and 1.50 r e ). Enough information is now available to 
surmise empirical radial profiles of the macroscopic parameters of such 
high-velocity streams. Figures 1, 2, and 3 show, respectively, typical 
radial dependences of the interplanetary magnetic field, the plasma density 
and the electron and proton temperatures. Several features of these profiles 
merit special emphasis: 1) the radial component of the magnetic field, B, 

between 1.03 r 0 and 2,r® decreases much more rapidly than the (l/r^) field 
of the spherical model; 2) the density above Coronal Holes is substantially 
lower than the average and exhibits a much steeper gradient; 3) the tempera- 
ture of the core electrons can be approximated roughly by an empirical poly- 
tropic relation: T c « n^ where 6 = .185 and n is the density; 4) the proton 
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temperature displays a strong correlation with the wind velocity (in par- 
ticular, the high-velocity streams are characterized by high temperature 
protons). This last point deserves further elaboration: Helios data show 
that the proton temperature, T p , in a high-velocity stream grows rapidly 
as one approaches the Sun and reaches about one million degrees at 0.3 A.U. 
This type of radial dependence of T p precludes any possibility of a poly- 
tropic relation between T p and n. 

The above sketched features of the solar wind beg for a physical ex- 
planation. As a first step in this direction, we have initiated a "feasi- 
bility" study of a model that emphasizes the role of the suprathermal elec- 
trons in the heat conduction along any given field line. We have attempted 
to adhere to first principles and the empirical facts as much as possible, 
and purposely ignored the possible importance of any other ad hoc driving 
mechanism. In other words, we have attempted to answer the following ques- 
tion: Do the suprathermal electrons generate enough heat flux to drive a 

high-velocity solar wind? The following is the preliminary report of our 
findings to date. 

2. Basic Equations for the Field-aligned Solar Wind Flow 

We confine our discussion to a steady field-aligned flow along a radial 
tube of force emanating from the central portion of a given Coronal Hole 
and concern ourselves only with the inner region of the solar cavity (from 
1.03 r Q to 1 A.U.). This allows us to ignore the effects of the rotation 
of the Sun and the curvature of the field lines. Furthermore, we assume 
that the r-dependence of B is known on empirical grounds. Finally, except 
for the appearance of a new term for the heat flux, we take the basic MHD 
equations as valid. Using standard notation, (B for the magnitude of the 
magnetic field, V for the wind speed, p for the mass density, P for the 
thermal pressure of electrons and positive ions combined, and q for the 
heat flux vector) we have for the plasma as a whole: 

pV = aB ( 1 ) 
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where a is a streamline constant. Eq. (1) follows from conservation of 
mass (div(pV)=0) and magnetic flux (div £=0). Equations (2) and (3) express 
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conservation of momentum and energy^ respectively.- Note that, since q is 
field aligned, one can replace div q by B*v (q-$/lr) or, for our geometry, 

div 5 = Hr pV ( 4 ) 

where 
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= _9_ 
ctB 


= _ 9 _ 
pV 


(5) 


If the function W is assumed to be a known function of the radial distance, 
r, or alternatively, of any of the dependent dynamical variables p, V, P, 
equations (1), (2) and (3) represent a closed system of equations for three 
unknowns: p, V and P. Two of these unknowns can be eliminated immediately. 

First, Equations, (2) and (3) yield the integral: 
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where U is a streamline constant representing the asymptotic value of the 
wind speed. Equation (6) allows us to eliminate P in terms of remaining 
variables. Secondly, Equation (1) allows the elimination of, say, V in 
favor of p. We obtain after some algebra, the following first-order differ- 
ential equation for p: 



r d p _ H(p ,r) 
p dr L(p,r) 

(7) 

where 
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L = 4(— ) 2 + 2W + 2p ^ - U 2 
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In Section 4, we shall explain the procedures which we used to con- 
struct an explicit form of W as a function of p. For the present we need 
only the result. Figure 4 shows the (dimensionless) quantity W/U 2 plotted 
versus electron density n ; it has been computed under the assumption that 
a 2 = 2.4 X 10-23 g/cm3, u e = 700 km/sec and that B of Figure 1 may be 
approximated by the empirical formula: 
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B = 2.2 (-p-)^/l+47 (^) 6 (gauss) 


( 10 ) 


where r = 1 .03 r . 
s ® 

For the purpose of numerical integration of Equation (7), we found it 
convenient to represent W in Figure 4 by an analytical formula: 



22 

where Pq * 10" g/cm 3 and W Q *10 4 (km/sec) 2 


3- Numerical Calculations of (p,V,P) profiles 

Detailed analytical study near the critical point of Equation (7) (i.e., 
the radial distance r c at which H=L=0) has shown that Equation (7) allows a 
two-branch solution through r c : one for the accretion and one for the 
outflow. The topology of the solution in the neighborhood of the critical 
point is analogous to that of the polytropic Parker model. Figure 5 shows 
the density profile corresponding to the outflow with U=700 km/sec. A 
comparison of Figure 2 data referring to the high-velocity streams (lower 
segments of the n e curves) reveals an excellent agreement with the curve of 
Figure 5. Once the density profile is obtained, it is a simple matter to 
compute the remaining profiles for V and P. Figure 6 shows the wind-speed 
profile based on Equations (1) and (10) and Figure 7 shows the quantity: 
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Since P - k(n T + n T + n T ) (n and T indicate the number density and 
the temperature of trie a-pa?t?cles, respe8ti vely) - T e ff may be interpreted 
as the effective" temperature of a high-velocity stream. A comparison of 
the T e ff -curve with the combined data for T c , T D and T a = 4T in the ex- 
plored region of the heliosphere (r > 0.3 A.U.) shows a surprisingly good 
agreement both in slope and the levels of the ion temperatures. The agree- 
ment is improved further when a (relatively minor) correction for the 
presence of the Alfven-wave pressure is made in the expression for P. Of 
course, the unexpected result of our calculations lies in the unexplored 

?Tin6 < r < T^f exhibits a rise with a maximum of 

Z.5xl0 K at 9 Since the core-electron temperature in our model de- 
creases monotonically with r (see next section) and thus at 9 r Q is sub- 
stantially lower than 2.5x10° °K, one is forced to conclude that the model 
requires surprisingly high proton and a-particle temperatures in the outer 
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regions of the corona. We must stress that if one accepts the form of the 
heat flux used in the calculations, this conclusion is forced upon us oy 
nothing else but the basic laws of conservation of mass, momentum and 
energy. Since our form of the heat flux predicts correct density profiles 
of the high-velocity streams, we find the conclusion inescapable. 

4. Heat Flux Calculations 


We now outline briefly the method and the steps that we have taken to 
arrive at the result shown in Figure 4. Space limitation does not allow 
us to go into detail but only to mention that one can develop a variety of 
cogent arguments demonstrating that the conventional theory for heat con- 
duction in collision dominated plasma (e.g., Spitzer's theory) is not 
applicable to any portion of the solar corona. The traditional Fourier law 
(q„ =-Kv„T) breaks down right at the base of the corona where T=10^ °K and 
n e = 108 cm" 3. This i s so i n S pite of the fact that the coulomb mean free 
path (mfp) of thermal _electrons, X c> at the base is still very small compared 
to the solar radius {x c ~ 10"^ r 0 ). The main reasons for it are: a) the mfp 
grows rapidly with the electron energy (a suprathermal electron of 7kT has 
a mfp 50 times larger than x c ); b) the steep density gradient above Coronal 
Holes leads to a short scale height H ( H - 0.1 r Q ); c) small H leads to a 
rapid development of skewness in the distribution function of suprathermal s; 
and d) the forward skewness of suprathermals is fostered by the paucity of 
collisions and the open line morphology of the rapidly decreasing magnetic 
field. For these reasons, any perturbation method relying on the smallness 
of the Knudsen number (x c /H) will lead to erroneous results concerning the 
physical characteristics of the heat flux. The reader is invited to read 
more on this in the paper presented by J. Scudder during these proceedings. 

There seems to emerge from these types of considerations a clear cut 
possibility: the suprathermals may carry enough heat to drive the solar 

v/ind. To check this in a quantitative manner we saw no way out but to 
rely on the first principles of the plasma kinetic theory, i.e., to evalu- 
ate first the electron distribution function, f e , at various points along 
some given tube of force and then by direct numerical integration compute 
the heat flux at these points. Since the solution of the Boltzmann 
equation for f e in its more rigorous form (say, in the Fokker-PJanck ap- 
proximation) is prohibitively difficult, we have opted for the Krook's 
approximation which, although not as accurate as the Fokker-Pl anck approxi- 
mation, is on physical grounds far superior to any perturbation method 
(such as that of Chapman and Enskog). 

First, let us list the simplifications that can be made in the 
general Boltzmann equation. We are interested in the distribution function, 
f e , of transthermal and suprathermal electrons. For these electrons, the 
proper plasma frame and the inertial heliocentric frame of reference 
practically coincide. The gravitational effects are negligible. Moreover, 
f e is gyrotropic, i.e., in the steady state, f e is a function of only three 
variables: radial distance r, the speed v and pitch angle e. With these 

stipulations, f e (r,v,e) obeys the following equation: 
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( 13 ) 
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where E„ is the component of the electric field parallel to $ and is re- 
lated to "interplanetary polarization potential", <j>, by E„ = -d$/dr; x is 
the mean coulomb collision time for momentum transfer and is given by 
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In our case the numerical factor lnA=25 and £ is of order of unity. For 
the "source function" f^, which coincides approximately with that of the 
distribution of the core electrons, we adopt a maxwellian form 


f M ( r ’ y ) 


n 

~VT 3 
r w_ 



(15) 


where w c is the thermal speed related to the temperature, T c , of the core 
electrons by n^w^ = 2kT c . Note that the number density, n e , and w are 
functions of r. Other forms for f]^, compatible with observations of thermal 
and subthermal electrons, may be chosen to improve the description of the 
electron transport. The reader will find more on this in the aforementioned 
report by J. Scudder. Since the brunt of the dynamical processes shaping 
the evolution of the solar wind takes place in the coronal region, we are 
primarily concerned about the adequacy of Equation (15) in that region. 

Since in the limit of x being very small f e should approach a maxwellian 
form. Equation (15) ought to represent a reasonable first step in an 
iteration process to which Equation (13) could be subjected. Note that in 
the limit of a very large x (collisionless plasma) Equation (13) has a well 
known general solution 


f e = *(E,w) 


( 16 ) 


where 4* is any appropriately chosen function of the two constants of 
motion: 
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(17) 




- e<f> » 


u = 


2 . 2 
nigV^ sin 0 

2B 


Not much is known about f e for finite t's« But it turns out that Equation 
(13) can be integrated formally for any given form of f^ and t by methods 
of characteristics. If one confines oneself to the particle trajectories 
in the forward direction (0 < 9 < 90°) one finds: 
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In Equation (20) v" is given by Equation (19) upon replacing r* by r" 


and 
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Note that, for a fixed r, v and 9, Equation (21) imposes a lower limit 
on r' given implicitly by the equation: 
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where v m is given by Equation (19) upon replacing r' by r m ; r,j, is the point 
where e" = 90°, i.e., the mirror point of the trajectory. This implies 
that care must be taken while evaluating Equation (18) that the inequality 
chain: 


r < r < r' < r" < r (23) 

o - m - 


155 



be not broken; r 0 identifies the boundary of the onset of the outward heat 
flow. The precise value of r 0 cannot be ascertained, but it ought to fall 
somewhere between 1.03 r 9 and 1.30 r 9 ; fortunately, f e for r > 2 r 0 is not 
sensitive to r 0 . The explicit numerical evaluation of Equation (18) re- 
quires the knowledge of n e , T c , 4) and B as functions of r. We made the 
following choices: 


n e (r) = exp[13.25 - 5.50 In £- + 3.50 V(ln - 0.82) 2 + .79] (24) 


V r > = T 1 ^ w 


+ (r) = wkT ' ;(u = 1+1) (26) 

C p 


B is assumed to be given by Equation (10) 


and 6 = 0.185, T ] = 8xl0 4 °K, n ] = 2 cm 3 . 

In constructing formulas (24) and (25) we were guided by the observational 
data referring to the high-velocity streams. Equation (24) reproduces well 
the Coronal Hole and the high speed profiles of Figure 2. Equation (25) 
agrees with the dashed curve in Figure 3 and the empirical findings of 
Sittler and Scudder; its validity in the solar corona must remain at the 
present moment an open question although the numerical values it offers in 
the inner corona are not in contradiction with the observational estimates. 
Equation (26) is primarily based on theoretical considerations involving 
the generalized Ohm's law; for our purposes, it is in a sufficiently good 
agreement with the findings of other authors; but again, its validity in 
the inner corona must remain an open question. 

The entire subject matter of finding self-consistent functional re- 
lations between np,T c and <j> from the first principles is still unresolved. 

It is a subject of current investigations by Dr. Scudder and myself. 

"Armed" with the empirical formulas (24) through (26), we are in a 
position to evaluate f~ as a function of v and 0, for various values of r. 
Since the computer costs involving this project have proven to be very high, 
we had to confine ourselves to only six values of r (r/r s =1,2, 4.7, 16, 

63 and 208). Figures 8, 9, 10 and 11 show the forward portions of the con- 
tours of constant f e plotted in the velocity space for r/r s =2,4.7,16 and 
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and 63. Note that the coordinates are v„ and v x measured in terms of the 
local thermal speed, i.e., the horizontal coordinate is v cose/w c and the 
vertical coordinate is v sine/w c . The ratio of the values of f e of two 
neighboring contours is /To"; the outermost (lowest value) contour corre- 
sponds to f e = 8 x10“3T (sec/cnr)3. As one proceeds outward from the Sun one 
witnesses a gradual evolution of a highly skewed, forward "strahl", very 
much like the strahl observed onboard Helios spacecraft. There are various 
other features in the distribution that are of potential value from the 
plasma-theoretical point of view; e.g., the development of a col in the 
transthermal region of the velocity space. Such cols are sources of various 
plasma waves which in turn affect the thermal electrons. Clearly, the re- 
sults shown here open up a new area of future investigations. 

We have not computed the backward portions of the distribution function 
(90° < e < 180°) for the simple reason that we believe that f e cannot de- 
velop any significant backward "strahl" and, for the purposes of the heat 
flux calculations, may be approximated by a maxwellian form for all speeds. 

The knowledge of f e allows us to compute the heat flux q at various 
values of r. We have 


q(r) = it m g ffv 


5 


cose f e (r,v,e) sine de dv 


(27) 


Because of the complexity of f e as a function of v and 0. the double 
integration required to obtain q is enormously expensive and we had to con- 
fine ourselves to the six r- values quoted above. Thus, the entire effort 
leaves us with only six numbers which must serve us as "beacons" for the 
construction of the curve shown in Figure 4. Needless to say, we had to 
caxe some iiDerties to periorm tne necessary inLer|ju1duiuri leauiny Lu tlic 
results shown in Figure 4. 

Finally, we should like to comment why we have chosen to represent the 
function W=q/(aB) as a function of density rather than as a function of 
distance. Firstly, purely on practical grounds, it is more economical to 
evaluate f e when one replaces the integration variable r' in Equation (18) 
by n e (r') or T e (r') via Equation (24) or (25). Secondly, the dependence 
of W on r is primarily predicated by its dependence on ne- Thus, on 
physical grounds, it is more instructive to know W as a function of n e 
rather than as a function of the arclength along a specific field line. This 
does not mean to imply that W may be looked upon as a unique function n e for 
all field lines. On the contrary, our calculations show that q is quite 
sensitive to the behavior of the magnetic field. Since the possibility for 
relating B to n e by some theoretical "law" is quite remote, we must treat B 
and n e as separate entities and thus regard W as a function of n e and B. 

5. Summary 

Our preliminary calculations show that the heat flux generated by the 
suprathermal electrons plays an important role in the dynamics of the solar 
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wind. In fact, there is a possibility that it might account entirely for 
the salient features of steady, pure high-velocity streams. Of course, 
more work is required to resolve the issue, but we believe that the initial 
results presented here invite such work as worthwhile. 


FIGURE CAPTIONS 

Figure 1: Typical radial profile of the interplanetary magnetic field, 

Except for th^ shaded area near the Sun, the curves are based 
on Parker's spiral model and are in agreement with observations 
averaged over fluctuations. 


Figure 2: Typical radial profiles of the electron number density, n e , in 

the solar wind. The curves were constructed from a variety of 
observational sources. The two lower segments labeled "density 
above polar Coronal Hole" and "fast streams" are characteristic 
of the central portion of a high-velocity stream. 


Figure 3: Typical radial profiles of proton and electron temperatures in 

the solar wind. The solid curves represent the scalar proton 
temperature, Tp=(l/3)(2T i + T„), for three adjacent intervals of 
the wind speed, I, II and III. The dashed curve refers to the 
core electron temperature, T_; the dot-dashed curve, labeled Th» 
represents 7T C and is meant to indicate the domain of the 
suprathermals. Except for the tentative markings near the Sun, 
the plots are a "smoothed out" version of the Helios and Voyager 
data . 


Figure 4: The heat conduction W, measured in terms U 2 and plotted versus 

electron number density, n e ; W is related to the heat flux q by: 
W = q/(pV) = q/(aB). See text for the definition of symbols and 
the theoretical considerations leading to this plot. (The arrow 
labeled C.P. indicates the position of the critical point). 


Figure 5: Calculated radial profile of the electron number density, n e . 

The curve represents the "critical" solution of Equation (7) 

(see text) and should be compared with the lower segments of the 
data of Figure 2. 


Figure 6: Calculated radial profile of the wind speed, V. The curve is 
computed from V = ctB/p using the result shown in Figure 5. The 
point labeled C.P. indicates the position of the critical point. 
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Figure 7: Calculated radial profile of the "effective temperature" of a 

high-velocity stream (see text for the definition of T e ff and P). 


Figure 8: Contours of constant electron distribution function, f e , plotted 

in the velocity space for r = 2r s . The horizontal (vertical) 
coordinate represents the parallel (perpendicular) component of 
the electron velocity measured in terms of the local thermal 
speed, vi~. The ratio of two neighboring contours is /TO; the 
outermost (lowest) contour represents f e = 8xl0" 3 ' (sec/cm 2 ) 3 . 
The figure implies gyrotropy and shows only the forward portion 
of f e (0° < 9 < 90°); r $ =1.03 r @ . 


Figures 9, 10 and 11: Same as Figure 8 but referring, respectively, to 

r/r $ = 4.7, 16 and 63. 
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